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1 Introduction 



Matrix models have been very useful in the study of the quantum geometry of two- 
dimensional quantum gravity. In [1] this program was extended to three-dimensional 
quantum gravity. It was shown how the so-called ABAB two-matrix model de- 
scribes the transfer matrix of three-dimensional quantum gravity.^ More precisely, 
a non-perturbative, background-independent definition of quantum gravity, which 
emphasizes the causal structure of space-time and which allows rotations between 
Lorentzian and Euclidean signature, was proposed in [2, 3], generalizing an explicitly 
solvable two-dimensional model with these features [4] . In the model, which has an 
UV lattice cut-off a which should be taken to zero in the continuum limit, one can 
define the concept of proper time. In the Euclidean sector the corresponding evolu- 
tion operator is defined in terms of the transfer matrix T describing the transition 
between quantum states at (proper) time n • a and (proper) time (n + l) ■ a. The 
transfer matrix is related to the quantum Hamiltonian of the system by 

f = e""^. (1) 
The ABAB model is defined by the two-matrix integral 

= ydyld5e-^*'-(^(^'+^')-^^'-'?^'-^^^^^). (2) 

Under the assumption discussed in [1] the free energy F{ai, a2, (3) is related to the 
matrix elements of the transfer matrix T in a way reviewed in the next section. 
The matrix model (2) has a scaling limit for ai = a-i which was analyzed in [7]. 
This allowed us in [1] to determine the corresponding phase diagram for the three- 
dimensional quantum gravity model and to map the bare coupling constants of the 
gravity model to the matrix model coupling constants a\ — OL2 and /? [9]. However, in 
order to study details of the scaling relevant to three-dimensional quantum gravity 
we have to study the matrix model for a\ ^ a2- In the scaling limit of interest 
for us both ai and a2 will scale to a critical value ac, but independently. Since we 
are interested only in the behaviour of the theory near the symmetric solution we 
need only the perturbative expansion around this solution rather than the complete 
solution in the asymmetric case^. 

The rest of this article is organized as follows. In sec. 2 we review shortly the non- 
perturbative definition of three-dimensional Lorentzian quantum gravity [2, 5] and 
its relation to the ABAB matrix model. In sec. 3 we review the machinery needed to 

^Previous work on 3D quantum gravity in terms of Lorentzian triangulations can be found in 
[5,6]. 

"^Whilc writing this article the asymmetric ABAB matrix model has been solved by Paul Zinn- 
Justin [10]. The behaviour close to the symmetric line ai = a.i is the same as the one reported 
here and to extract it one has to expand the elliptic functions which appear in the solution, an 
effort comparable to the one used here. 
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solve the ABAB matrix model for a symmetric choice of coupling constants [7]. In 
Sec. 4 we discuss the solution of the general ABAB matrix model, and in Sec. 5 we 
expand around the symmetric critical point relevant for three-dimensional quantum 
gravity. In Sect. 6 we discuss how to extract information about the transfer matrix 
of 3D gravity, knowing the free energy of the asymmetric ABAB matrix model. 

2 Quantum gravity and the ABAB matrix model 

Simplicial Lorentzian quantum gravity in three dimensions is defined in the fol- 
lowing way: the spatial hypersurfaces of constant proper-time are two-dimensional 
equilateral triangulations. Such triangulations define uniquely a two-dimensional 
geometry. It is known that this class of geometries describes correctly the quantum 
aspects of two-dimensional Euclidean gravity. It is also known that the description 
of two-dimensional Euclidean quantum gravity in terms of the class of (general- 
ized) triangulations is quite robust. In [2] we used this universality in the following 
way: the two-dimensional geometry of the spatial hypersurfaces is represented by 
quadrangulations and it was shown that it is possible to connect any such pair of 
quadrangulations by a set of three-dimensional generalized "simplices" . More pre- 
cisely, let a be the lattice spacing separating two neighbouring spatial hypersurfaces 
at (propcr)-times t and t + a. Then each square at t is connected to a vertex sXt + a 
and each square at t -|- a is connected to a vertex at proper-time t, forming pyra- 
mids and inverted pyramids. A further needed three-dimensional building block is a 
tetrahedron connecting a spatial link at t to a spatial link at t+a. The proper-time 
propagator for (regularized) three-dimensional quantum gravity between two spatial 
hypersurfaces separated by a proper time T = n-a is obtained by inserting n— 1 inter- 
mediate spatial hypersurfaces and summing over all possible geometries constructed 
as described above. The weight of each geometry is given by the Einstein action, 
here conveniently the Regge action for piecewise linear geometries. According to 
[3, 5], the contribution to the action from a single discrete time step is given by 

S = -K{Nt + Nt+a - N22) + X{Nt + Nt+a + \n22), (3) 

where Nf, Nt+a ^-nd N22 denote the number of pyramids, upside-down pyramids and 
of tetrahedra contained in the slice [t, t + a], and k and A are the dimensionless bare 
inverse gravitational and bare cosmological constant in three space-time dimensions. 
The naive continuum limit is obtained by scaling the lattice spacing a ^ while 
keeping T = n ■ a fixed. However, different scaling relations between T and a might 
in principle be possible^. 

*In two-dimensional Euclidean quantum gravity the proper-time T scales anomalously and one 
has to keep n^/a fixed. This is in contrast to the situation in two-dimensional Lorentzian quantum 
gravity as defined in [4] where the proper-time T scales canonically. The relation between the two 
models is well understood [11]. 
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Let gt and gt+a be spatial two-geometries at t and t+a, i.e. two quadrangulations 
and let {gt+a\T\gt) be the transition amplitude or proper time propagator from t to 
t+a. By definition, T is the transfer matrix in the sense of Euclidean lattice theory, 
and it satisfies the axioms of a transfer matrix [5]. In the case where the spatial 
topology is that of it was argued in [1] that the continuum hmit could be obtained 
as the large N scaling limit of the matrix model (2). Let Nt and Nt+a denote the 
number of squares in the quadrangulations associated with gt and gt+a- The two- 
volumes of the corresponding geometries are thus NfO^ and Nt+a(i^, respectively, and 
the relation to F{ai,a2, (3) defined by (2) is 

F{a,,a2,f3)= ^ e-^*^*(«*)-^*+''^*+'^(^*+'')(^i+„|f |^,), (4) 

gt,gt+a 

where Zt and Zt+a are dimensionless boundary cosmological constants. The naive 
relation between the matrix model coupling constants and the bare dimensionful 
gravitational and cosmological couphng constants and A^'^) of three-dimensio- 
nal gravity is 

ai^e^-^''*, a2^e''-^~'*+% /3 = 6-^5^+'^), (5) 

where 

In this paper we shall discuss the non-perturbative renormalization of the coupling 
constants. In principle we are interested in the limit Zt — Zf+a — 0, i.e. ai — a2- 
However, in order to be able to extract the information about the scaling of the 
boundary cosmological constants we have to keep Zt and Zt+a different from zero at 
intermediate steps. Thus these boundary cosmological constants should be viewed 
as source terms for the boundary area operator. 

3 The symmetric case: a.i = 0.2 = ct 

Let us for later convenience shortly review the technique for solving the matrix 
model (2) used in [7] (based on earlier results [8]). 

By a character expansion of the term g^^^*''^^"^^ one can write 

Z{ai,a2,P) ~ ^ (^)' ^ R{h}{ai)R{h}{(^2), (7) 

{h} 

where the sum is over the representations of GL{M), characterized by the shifted 
highest weights hi = nii + M — i, (i = 1, . . . , M), where the are the standard 
highest weights and where the large-M limit of the coefficient c^h} is 

log^M = - E|(l°g^ - l) - ^log^(^)' ^(h) = Hih. - h,). (8) 

i i<j 
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Finally, if X{h} denotes the character associated with {h}, 

%}(«) = / MhM) expM( - ^tvA^ + jtrA'^y (9) 

It is now possible to perform a double saddle point expansion of (7) and (9). In 
order to describe the formalism let us introduce the notation 

^^(^),/(>~ + «) + /(^-.o)_ a/w = ^'^ + '°);^''-'°'. (10) 

This notation is useful when f{z) has cuts. The saddle point expansion assumes the 
existence of an eigenvalue density p(A), or equivalently a resolvent associated with 
the matrix integral (9): 

^(^) = 7^Ea3Y' -^im-^^w, (11) 

k 

and (after rcscaling h — > h/M) a density of highest weights p{h), or the correspond- 
ing "resolvent" H{h): 

^(^) = ]^E/^' -7Ttp{h) = QH{h). (12) 

k 

In [7] the double saddle point expansion is analyzed in the case ai = a2 = a. 
The density p{h) was assumed to be different from zero in the interval [0,/i2[, and 
equal to 1 in the interval [0, hi], where < hi < /i2- Further, for a given eigenvalue 
distribution of the matrix A coming from the saddle point of (9) one can define 
a function L{h), with same cut as i7 by 

mhj) ^ ^-^^logXiH}{A{X,)). (13) 

The analysis of [7] shows that L{h) = H{h) + F{h) where F{h) is analytic on the 
cut of H[h) but has an additional cut [/la, oo[ where 

2^L(h)^ log- + H(h). (14) 

a 

It can now be shown that the function D{h) = 2L{h) — H{h) — 3 log h + \og{h — hi) 
only has square root type cuts on [hi, /i2] and [h3,oo[ and on these cuts satisfies the 
following equations: 



^D{h) = log^^;^' heIo^[hi,h2] (15) 



m{h) = log^^, heIi^[h^,oo[. (16) 
5 



-i -1 i +1 
1 

hi +1 hg -1 

Figure 1: The cut structure of r(/i) in the complex /i-plane. 

Eqs. (15)-(16) constitute a standard Hilbert problem and the inversion formula is 
unique [12] . The function holomorphic in the plane with cuts Iq and Ii is given by 

D{h) = log-— i---°^r(/i) / dh'j———— (17) 



Ph? in ' 'A3 {h-h')r{h') 



+^(^) / d/i'-- -2r(/t) / d/i' 



J —ex. 



{h-h')r{h') ' 



where 



r{h) = ,J{h-h,){h-h2){h-h^) (18) 

and where we have chosen the cut structure shown in Fig. 1. Following [12] the 
meaning of r{h') on the cut is r(/i'+iO), i.e. the function on the "left side" of the cut. 
The integrals can be expressed in terms of standard elliptic functions. However, we 
do not need the exphcit expressions here. 

Prom D{h) we can derive the expression for p{h) which is 

p{h) = = h e]hi, h2[. (19) 

ITT m 

We have {h always h+iO if ambiguities) : 

-r{h) r^' dh' _ {-2r{h)) dh' 

~ ~^J-oo {h-h')tr{h') n {h - h')tr{h') ^ 

log/3/a, ,,,, dh' , , 

Note that the derivative of D{K) and p{K) with respect to hi are elementary 
functions of h. For instance, differentiating with respect to /la we have 

dD{}i) _ ir{h) dW _ niJi)F:, _ . dp{h) 
dhs ~ h-hs dh^ ~ 2{h - /la) ~ dh^ ' ^ ' 

where the last equality is valid for h E Iq — [/ii,/i2] and where W{hi,h2,hs) and 
Fz{hi,h2,h'i) are defined below (eqs. (25) and (28)). Thus we can write 

D{h-M-r5h:,) = D{h;hs)+Shs'-^^^^^ (22) 
and similarly for p{h). The function F^lhi, /i2, h^) is a sum of elliptic integrals. 
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3.1 Boundary conditions 



The starting formula for D{h) is (17). The general large-/i behaviour of this function 
is 

cih^/^ - log{-ah) + C2/i-^/2 + 0(1//?,). (23) 

However, according to the analysis in [7] Ci = and C2 = — (— a)^^/^. This gives two 
boundary conditions for the constants hi,h2,h3 which appear in r{h) and thus in 
D{h). The coefficients ci and C2 can be identified by expanding the integrand in 
(17) in powers oil/h and one obtains the boundary conditions: 



ci = iW{hi, /i2, ^3) = 0, C2 — i^ — 




where we have used the first of the equations (24) to simplify the second, and where 
,,,,, , , , \og(5 la r d/i' dh' dh' , , 

, , , log (5/ a f^^ h'dh' h'dh' ^ f^' h'dh' 

We will not need the explicit expressions for these integrals. The final boundary 
condition is 

/ dh p{h) = 1 - /ii, (27) 

where p{h) is given by (20). 

For given (a, /5), we can in principle solve the three boundary conditions (24) and 
(27) for (/ii, /i2, /is), thereby obtaining a solution of the matrix integral. Equivalently, 
the three conditions define locally a map from the /3-Q;-plane to a two-dimensional 
hypersurface in the parameter space of the h^. Critical regions of the free energy F 
of the model are associated with singularities of the inverse map of an independent 
subset of the /ij's, say, (/i2,/i3) ^ {o:{h2, h^), P{h2, hs)), which lead to singularities 
of F{a,l3) upon substitution. 



3.2 The critical line 

The generic behaviour of D(h) when h is clearly (h — h^)^^'^, simply coming 

from the term r{h) in the representation (17). However, this behaviour can change 
to {h — hsY^'^ along a curve adP) in the (/3, a) coupling constant plane. According 
to [7] this is the critical line of phase A of the ABAB matrix model and according 
to [1] this is where the continuum limit of 3d gravity should be found. Similarly 
the criticality in the B phase is derived from the behaviour of D{h) when h — > /i2, 
where a generic behaviour {h — h2Y^'^ changes to {h — /i2)^^^- 

We now study the change of (/ii, /i2, h^) as a and (3 change infinitesimally. For 
simplicity we first present the result when a/ (5 is constant. 
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Let us first identify tlie coefficient of \/h — in D{h). Using (21) in the expres- 
sion for D{h) tlie coefficient can be written as 

iro{h3)F3{hi,h2,h), ro{h) = - - /la), (28) 

One lias 

^55%^ = iF3(ft.,,>../,3). (29) 

We define Fi and F2 similarly to F3 and have relations like (29). From the boundary 
conditions (24) it follows that the variation of hi, h2, as a, (3 change with the ratio 
a/ j3 fixed satisfy 

Fi6hi + F25h2 + Fsdh = 0. (30) 

hFi6hi+h2F26h2+hsFs6hs = -2^. (31) 

The final boundary condition involves the density. Since p{hi) — l and p{h2) — 
the variation of (27) just becomes 



+ + (32) 

Jh^ \dhi dh2 dhs ) 



which can be written, using (21), as 

EiFibhx + E2F2bh2 + E^F^bh^ = 0, (33) 

where 

.^1.2.3 (34) 

fel tli — tl 

are elliptic integrals. It is straightforward to repeat the derivation in the case where 
the ratio a/ (3 is not assumed constant, leading to the set of equations 



+ hiFi5hi+h2F25h2+h3Fs5h3 = 0, 



Sa 


Sf3 


a 


~J. 


Sa 


Sf3 


a 


J 


Sa 


S(5 


a 





a3/2 

R3 [^-^)+ EiFiShi+E2F2Sh2+EsFsShs = 0, (35) 

for the linear variations, where the functions Ri{hi, h2, h^) are given by 

2 dh' 2 /■'^^ h'dh' 

2 /•°° (]],' 



Without loss of generality, we can now choose /i2 and as independent variables 
on the 2d hypersurface. Using (35) to eliminate 6hi, one can derive the associated 
linear map between the remaining variables, 




X1F2 Y1F3 

X2F2 ^2-^3 




Sh2 



(37) 



where Xi and Yi are easily calculable functions of the i?'s, £"s and /I's. This trans- 
formation is in general well defined, unless either F3 = or F2 = 0, making the 
Jacobian vanish. We already know that these two equations define a critical line in 
coupling- const ant space, the former giving rise to phase A, and the latter to phase 
B. They meet in a single critical point with ^3 = ^2 = 0. 

In determining a continuum limit of the matrix model related to three-dimensional 
quantum gravity, we are interested in the relation between (/9, a) and (/i2, ^3) along 
specific curves as they approach a point on the critical line. Note first that along 
generic curves and away from the critical line, by virtue of (37) all variations will 
be of the same order, namely 



This behaviour changes when the critical line is approached from an infinitesimal 
distance. As can be seen from (37), any such curve whose tangent does not coincide 
with that of the critical line in phase A at their mutual intersection point has 
vanishing derivatives da/dh^ and dj3/dh^ there, by continuity. This means that 
although the variations 5a and 5/3 are of the same order as 5h2 (and therefore also 
5/ii), their relation with 5h^ is of higher order, indeed. 



A completely analogous statement holds in phase B of the model with /i2 and 
interchanged, since the critical line is defined by F2 = there, instead of = 0. 
The qualitative difference between the two phases will only become apparent in the 
discussion of the asymmetric case below. 

An alternative way of approaching a point on the critical line that turns out 
to be relevant for quantum gravity is by moving along the line itself. This case 
is analyzed most transparently by adding (in phase A) the constraint -F3 = and 
expanding it along with (24) and (27). One verifies by a short calculation that the 
first-order variations in this case behave according to (38). 

4 The asymmetric case cxi ^ 0.2- 

As mentioned earlier, the construction of the transfer matrix requires that we perturb 
away from ai — a2- Let us discuss the general structure of the matrix model with 
7^ q;2 (as explained above we will only need an infinitesimal perturbation away 



5q; ~ (5/9 ~ Sh2 ~ Sh^. 



(38) 



(39) 
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from ai — a2 in the continuum limit) ^. The main difference in the analysis of the 
matrix model with ai 7^ a2 compared to the situation ai = a2 = a is that the saddle 
point solution involves two eigenvalue densities pi(A) and P2(A) corresponding to 
the two one-matrix integrals (9) with a — ai and a — a2- Similarly, we will have two 
functions Li{h) and L2{h) corresponding to (13) since the eigenvalue densities pi(A) 
and P2(A) appear via the matrix A{X) in (13). On other hand we have only one 
density p{h) coming from the saddle point of (7). In order to solve the saddle point 
equations it is natural to follow the same strategy as in [7] and make an educated 
guess about the analytic structure of the functions involved and then show the self- 
consistency of the solution. Since we have two functions Li{h), associated with the 
same p{h) but different p(A)'s, and the cut of L{h) from [/i3,oo[ can be traced to 
the saddle point equation for p(A) (see [7] for a discussion), it is natural to assume 
that Li{h) and L2{h) have a cut from [0, h2] (with p{h) = 1 in [0, hi]), and that they 
have separate cuts [hf\oo[ and [hf\oo[. In the case where (3 = this structure is 
indeed realized. 

We can now write down the generalization of (15)-(16) 

^Di{h)] = log^?^ + 3f?[/(/i)] yheIo = [hi,h2] (40) 



nDi{h)] = log^^ \/heh = [h^\oo[ (41) 



and 



^D2ih)] = log - 5R[/(/i)] yhelo = [h,h2] (42) 

^D2{h)] = log^^^ yhel2 = [hf,oo[. (43) 

In (40)-(43) the D's are related to the L's and the function H as below eq. (14), 
that is, 

Di{h) ^ 2Li{h) - H (h) -3 log h + log{h- hi), i = l,2, (44) 

where the subtractions of the log's are made to ensure that the functions DiS have 
square root cuts. As in the case of a single a, we assume that Li{h)—Fi{h) +H{h) 
where Fi{h) is analytic on the cut Iq of H{h).^ The function 

f{z) = Fi{z) - F2{z) = Li{z) - L2{z) (45) 

is at this point unknown, but we can write ^[f{z)] = f{z) on Iq. If we assume f{z) 
is known on Iq, eqs. (40)-(41) and (42)-(43) arc standard singular integral equations 

^As already mentioned an explicit solution of the asymmetric ABAB model has been published 
while this manuscript was being completed [10]. 

^The Fi should not be confused with the functions of the same name defined in (29) above. 



10 



of the Hilbert type and can readily be solved and one can write 



(46) 
(47) 



where D'l\{z) are given by formula (17) with (a, /is) = {ai,h^^) and (a, /is) 
{a2,hf^), respectively, and 

rm = ^J{h-h,){h-h2){h-hf). 

Furthermore, we have 

'^[Dk{h)]^-mp{h), A; = 1,2. 
Therefore, the "imaginary" parts of eqs. (46) and (47) are 

dt m 



mp{h) = nipWh)-n{h)-l 

Jio [t - h)ri{t) 

i7ip{h) = mpf{h)+r2{h)-l'^'^ ^^^^ 

Jin 



ho {t - h)r2(t) 



yh e lo 
V/i G Jo 



(48) 
(49) 

(50) 
(51) 



where f- is the principal value of the integral. Eqs. (50)-(51) determine f{z) and 
p{h) in terms of the densities p\^,p\^ (eq. (57)), corresponding to a — ai, h^ — h^^^ 
and a = a2 , h^^hf\ We can obtain an equation for f{z) by subtracting (50) and 
(51), leading to 



m{p\%h) - pl^h)) = i 



■ dtfit) fn{h) , r2{h) 



+ 



m{t - h) \ri{t) r2{t) 
4.1 Uniqueness of the solution 



(52) 



Let us discuss the solution of (52) which is a singular integral equation. In order to 
bring it into a standard form of singular integral equations, and for convenience of 
later applications, we introduce the notation 



^p'^{h)^'-^[p\\h)-pl\h)), 



(53) 



and 



{t - h) k{h, t) = - 



t-h ( 


h - 






h- 




h-h3y\ 


t - 




\ 


t - 





h-h3 
t-h. 



(54) 



11 



where from now on will always refer to the average 

h3 = l{h^^ + hf\ (55) 

and the function r{h) will refer to (18) with given by (55). 

The function k{h,t) is regular at h — t. Let us further introduce 

We can now write eq. (52) as 

/ - ^ + / - k{h,t)m = Ap^ih), (57) 
Ji^ m t-h Ji^ m 

where only the first integral is singular. The so-called dominant part of this singular 
integral equation is given by 

*M = A,,(,). (58) 

hi t — n 

This equation has precisely one zero mode, namely, 

^ ^/ ^ ^ JL (59, 

ro{t)Ji^m h-t ro{t) 



where ro(t) = \/ {t — hi){t — /i2)- Expressed in terms of fit) we have 

/<')^v^tt§(.^-)- 

By moving the /c-term in Eq. (57) to the right-hand side, we can repeat the steps 
leading to (60) , and moving the /c-term back to the left-hand side we finally obtain 

where the kernel N{t, s) is a Fredholm kernel, 



Vt^I^AHk^^ (62) 

In general the solution to the Predholm equation (61) will be unique [12]. We thus 
have a one-parameter family of solutions fc{t)- 

In order to determine the /our parameters hi,h2,h^\hf^ and the constant C 
we need four boundary conditions and one more condition, which in this case will 
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be the normalization condition for p (which appeared above in formulas (50) and 
(51)), 

hi+ p{t)dt -1^0. (63) 

Jhi 

The boundary conditions are again obtained by the requirement that the large-/i 
asymptotics of Dj{h) contain no h^^^ term while the coefficient of the h"^^^ term 
equal (— ccj)"^/^. This results in four boundary conditions 

Wr.= w,-i-iylf{m^o, (64) 

where i = l,2, Wj^W{h,h2,h^^^) and Qj^Q{hi,h2,h^f). 

The set of constraints (63) (64), (65) implies that in the asymmetric case the 
theory is defined on a three-dimensional hypersurface of the eight-dimensional pa- 
rameter space («i, a2, (3, C, hi, h2, h^3 \ hf^). Singular points of the map between the 
parameters {ai,a2,P) and some other parametrization of this hypersurface (say, in 
terms of (/t2, hf^), after eliminating C and hi) will correspond to critical points 
of the theory. 

In analogy with the symmetric case we expect criticality to be related to a change 



in the behaviour of the Dj{h) when h approaches one of the hi^^ (phase A) or h 
(phase B). In the first case, this is signalled by the vanishing of one of the coefficients 
jFg"'^ in the relevant expansion 
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Dj{h) = Ti^^ \Jh- hf + higher order, (66) 
with 

Ji,^{h^^-t)r,{t) 

where Fg"*^ = F^lh-^ = h^^^). The two critical hypersurfaces in phase A defined by 
\ j = 1, 2, intersect along a critical line which coincides with the critical line of 
the symmetric model. 

In phase B, there is a single smooth critical hypersurface defined by the simulta- 
neous vanishing of the two functions ^2^ (defined similar to (67)) appearing in the 
expansions 

Dj {h) = T^i^ \/h - h2 + higher order. (68) 

The simultaneous change in behaviour of t'^^ and T'^^ can be traced to the fact 
that the non-analyticity at h2 is due entirely to H(K) (or p(/?.)) which arc common 
to Di{h) and D2{h). The critical hypersurface defined by ^^2"'-' = intersects both 
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hypersurfaces J^^^ = and ^3^^ =0 in critical lines, and there is a single (tricritical) 
point where all of the three hypersurfaces meet. 

One can now proceed as in the symmetric case by expanding the five constraint 
equations to first (or any desired) order, and eliminating the dependent variations 
5hi and 5(7, say. The three remaining first-order equations can be solved for 5ai, 
Sa2 and 5(3 and put in matrix form. 



/ 6ai/ai\ 

5a2/a2 
V 5(3/ P J 




\^a2/a2\= X2iTf^ + Z2i X22J'f^ + Z22 X2^T't^ + Z2^ \5hf\. (69) 

where the functions X^^ and Z^^ are regular and generically non- vanishing. Moreover, 
it can be shown that the three functions Z^x vanish at points h — h^^ of the critical 
surface defined by J^^^^ — 0, by virtue of the fact that the derivative df/dhf^ at such 
points is proportional to the zero-mode of /. (Analogous statements hold for Zi2 
and Zi^.) The situation is therefore very similar to what happened in the symmetric 
case, cf. (37). Namely, the rank of the Jacobian of the coordinate transformation 
(cti, a2, (3) H- > {h^^\ hf \ /i2) drops by one when we go to one of the critical surfaces. 

Expressed in terms of the linear variations, this means that away from any of 
the critical surfaces one has 

5ai ~ 5a2 ^ 6(3 ^ 6h2 ~ Shf^ ~ 6hf\ (70) 

and that this behaviour changes when a critical surface, say, J^^^^ = is approached, 
to 

Sai ~ Sa2 ^5(3^ Sh2 ~ {Sh'i^f ~ Shfl (71) 

Similarly, it follows that if the critical (symmetric) line of phase A is approached 
along any path that is not tangential to either of the critical surfaces j — 1,2, 
the variations will behave like 

5ai ~ 5a2 ^5(3^ 5h2 ~ {Shf^f ~ {Shff. (72) 

As before in the symmetric case, if one chooses a linear variation tangential to 
one of the critical hypersurfaces, say, ^3^^ = 0, this imposes a linear relation between 
5h^^\ 5hf^ and 5h2, and at the same time forces 5h^^^ to be of the same order as 
the other two, leading again to the behaviour (70). The same is true if we consider 
a variation tangent to the critical line in phase A. This case is analyzed easily by 
adjoining the two additional conditions J-'^''^ = to the other boundary conditions, 
and expanding them to first order. 



5 Expanding around cxi = CX2 

In applying the matrix model to 3d quantum gravity, we do not need the complete 
explicit solution of the asymmetric ABAB model, but only infinitesimal variations 
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around symmetric points with a\ = and h^^^ = h^^^ . In this context, it is convenient 
to work with the symmetrized and anti-symmetrized variables 



a^y/^^, a=J — , /i3 = ^2_^-^, /i3 = ^3-_^. (73) 




Infinitesimal linear variations around any given point (a(o),a(o),/3(o)) or the corre- 
sponding {hf\ hf \ h2^) will be denoted by 6a, Sa, 6(3, 5h^, Sha and 6h2. In what 
follows, we will study further the relations between these variations which follow 
from the conditions (63), (64) and (65), since they determine the continuum physics 
of the model. 



5.1 Finding f{t) and Di^2{h) 

The discussion in this section will concentrate on phase A, which is the most relevant 
one for the quantum-gravitational application. We will comment briefly at the end 
on what happens in phase B. 

In order to make maximal use of the relations we already derived for the sym- 
metric case, it is convenient to perform a partial expansion of the relevant functions 
around points in the symmetric plane (characterized by a = 1, or, equivalently, 
113 = 0). Wc will Taylor-expand around hf^ = and simultaneously allow for a flnite 
shift to a = e^ away from the symmetric value a*-"-' = 1. 

Recall now the integral equation (61) which determines f{t). Under the variation 
introduced in the last paragraph, because of symmetry in h^^^ and h^^\ the function 
k{t,h) has an expansion in even powers oi dh^, 

k{t, h) = {Shfh{t, h) + {Sh3)%{t, h) + ---. (74) 

Thus knowing Ap''^{h) allows us to calculate f{t) perturbatively in Shs. We can 
also expand Ap^^{h) in Sh^ and e. This makes the integration on the left-hand side 
of eq. (61) possible order by order in terms of elementary functions. The expansion 
of Ap^^{h) is based on the following two observations. Firstly, we have for any (not 
just infinitesimal) e the exact relation 

D'''{h; ae') = D'''{h; a) + e g{h), (75) 



. LrW ^^-1^. (76, 

Next, D^^ih-jh^) and Q{h) have expansions of the form (22) with respect to their 
/i3-dependence. Thus l\p^^{h) can be written as 

Ap^^{h; Shs, e) = (oi Shs + 61 e) + (02 Shs + 62 s) {Shf H , (77) 
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where only 61 is not an elementary function of h. Explicitly, we have 

^p'%h) = (^-Sh - e g{h, h)^ +■■■. (78) 

Here and in the following we use the shorthand notation F^lhs) for the function 
Fslhi, h2, /i3, a, P). The integral in (61) can now be performed and one obtains 

f{t) = (c'^t^, + 5hs 4l^F3(/i3) -e)+---, (79) 



where C = C + 5hsF3{h3,) + J^^ -^j^- Note that until this point we have not 

made use of any of the boundary conditions, and therefore could treat 5h^ and e as 
independent. 

Finally we can insert f{t) in (46) and (47) to obtain Di^2{h). To order {5h^Y we 
obtain 

Dj{h; hi, h2, hi; aj, (5) = D^\h; h^, /i2, h^; a, f3)+ (80) 

-r iro{h^)F^{hri) 



± C [yh - /is + 0(5/13)) ± 5/i3 



2s/h - hs 

{5hsf iro{h) fdFM , Fs{hs) , F,{h) , FM 



4 Vh^HTs V dh 2{h - hs) 2(/i3 - hi) 2(/i3 - /is) ' 
Note first that C' — O since the large-/i asymptotics 

Dj{h) ~ mjh^/^ - logajh + iOjh-^/^ + 0{l/h), (81) 

must be satisfied for both D^'s to zeroth (and, of course, to any) order in Sh^. 
Furthermore, note that the ^-dependence to this order is simply due to a shift in 
the argument from logo; to logctj = logci; ± £ in the terms log((/i — hi)/{h'^aj)) in 

In order to give a more detailed discussion of the continuum limit, we will also 
need the explicit form of the boundary conditions (63), (64) and (65) to second order 
(5^3)^. From (49) we get 

^ %Di{h) + D2{h)] = -mp{h), h e lo, (82) 

which implies that 

p{h) = p^'ih; hi, h2, hs; a, (5) + 0{{5h^f), (83) 

where the explicit form of the 0((5/i3)^)-corrections can be read off from (80). An 
expansion of the constraints VVj to second order yields 

Wj = W^\hi,h2,h^,a,(5) (84) 
{5h,Y f dFM Fsih) Fsjhs) \ ^ 
^4 \ dhs ^2(/i3-/ii) 2(/i3-/i2)V 
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In a similar way we obtain explicit expressions for Oj to order {ShsY. 

These expressions can now be used to refine the relations (72) which govern 
the behaviour along an arbitrary, non-tangential approach to the critical line. One 
obtains 

So r^Sp ~ 5h2, 

5hz ■ Shs ~ Sa, (85) 

for some real a^. These relations can be diagonahzed to give 

{Shs ± Shf - {Sh'i^y ~ aiSa + os^/S ± aaSa, (86) 

characterizing the critical behaviour related to the approach to the two critical 
surfaces discussed earlier. 

By contrast, in phase B we have only a single smooth critical hypersurface charac- 
terized by ^g"*^ = 0. Again, moving tangentially to this surface will induce linear vari- 
ations which are all of the same order, whereas generic, non-tangential approaches 
to the surface will be characterized by 

6a 6(3 ^ da ^ 6h^ ~ 6hy, ~ (5/i2)^- (87) 



5.2 The scaling relations and renormalization of 3D gravity 

We now want to relate the scaling limit of the matrix model discussed above to 
the continuum limit of the discretized 3D quantum gravity. (Euclidean) quantum 
field theories can be defined as the scahng limit of suitable discretized statistical 
theories. The continuum coupling constants are then defined by a specific approach 
to a critical point of the statistical theory. Different approaches to the critical point 
might lead to different coupling constants or even different continuum theories. In 
our case we want to show that it is possible to approach a critical point in such away 
that the canonical scaling expected from a theory of 3D gravity is reproduced. 

We saw above that the asymmetric model is defined on a three-dimensional hy- 
persurface in the parameter space of (/ii, ^3^)- Let us consider a curve approach- 
ing a specific symmetric point {h\^ h"^^ h^^^ = ^3) on the critical surface of the model, 
with curve parameter a (where a = corresponds to the critical point). In terms 

of independent parameters, this path can be described by (/i2(a), h^i\o)i hf\ci)) or, 
equivalently, {ai{a),a2{a), P{a)). Recalling our discussion in Sec. 2 above, we are 
interested in obtaining a behaviour of the form 

log ai = log etc + ^TT" ~ ^^^i ~ a^C2A (88) 

Qj C3 

log/3 = log /3c - -pr - a^C4A 



17 



as a — > 0. In (88), logttc and log/?^ represent combined additive renormalizations of 
all the coupling constants, and G, A and Zi are the renormalized gravitational and 
bulk and surface cosmological coupling constants. 

For the purposes of conventional quantum gravity we are primarily interested in 
phase A of the model. Moreover, because of the symmetry between in- and out-states 
which is related to our three-dimensional geometric interpretation of the essentially 
two-dimensional matrix model, we are only interested in continuum limits where 
both of the discrete spatial boundary volumes go to infinity, that is, where both zt 
and ZtJ^a are critical. Interesting critical points therefore must lie on the critical line 
where the two critical surfaces intersect. 

As we have argued at length in [9], a standard critical behaviour of 3D gravity 
governed by the canonical scaling of the three-volume can only be achieved by a 
tangential approach to the critical line (in that case, in the ^-a plane). This ap- 
proach ensures that the terms proportional to a/G in the expansions of a. and (5 by 
themselves do not determine the leading critical behaviour, although they contribute 
in the combination AG^ at higher order. The same construction carries over to the 
asymmetric case, which differs from the symmetric situation by the appearance of 
the second-order terms proportional to Zj in (88). As happened there, the tangent 
vector to the critical curve (in the hyperplane a = 0) at a given point q;'^(/3) fixes 
the ratio C1/C3 of the constants appearing in (88). Moreover, the leading critical 
behaviour will now be governed by the boundary cosmological terms. This means in 
practice that in order to determine the singular behaviour of the free energy Ficii. /?) 
under such an approach, we can simply ignore the order-a terms in (88) and use the 
scaling relations (85) for the remainder. 

We still have to discuss the relative scaling behaviour oiba. and 8 a. As we will 
show in the next section, if we want the transfer matrix to reduce to the unit operator 
in the limit a ^ 0, the singular behaviour should to leading order depend only on 
ba and not on ^a. One way of realizing this possibility would be by showing that 
at this order the coefficient in front of the ^a-term vanishes, implying the relations 

ba oc a^, (89) 
boL oc 

and consequently 

5/i3 OC a, (90) 
(5/i3 oc c? . 

The corresponding scaling curve deviates from the symmetric plane a\ — a2 only by 
0(a^)-terms. Note that this means that the asymmetry 6a contributes at the same 
order as the cosmological term A, which presumably is a desirable property in view 
of the construction of the Hamiltonian. 
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6 The transfer matrix 



The transfer matrix contains the information necessary to derive the Hamiltonian 
of 3D gravity. Prom the free energy of the asymmetric ABAB model we can extract 
some information about the transfer matrix as is clear from formula (4). Let us 
be more precise about this (see [1] for a detailed discussion). The free energy of 
the asymmetric ABAB matrix model involves according to (4) a summation over 
the individual geometric states \g) which label in- and out-states. However, one 
can use the free energy to extract information about the areas Nin and Nout (the 
number of squares in the in- and out-quadrangulations) of the in- and out-states 
\gin) and \gout)- We expect this quantity to capture the essential part of physical 
information about the time evolution of a two-dimensional universe (cf. e.g. [14]). 
Let us consider the state 

l^) = ;^^^^EW*)l^t> (91) 

where J\f{N) is the number of quadrangulations of given area N. The norm of such 
a state is 

{N'\N) = dr,,N' (92) 

since states \gi) and \g2) with different quadrangulations are orthogonal. The num- 
ber of quadrangulations constructed from squares grows exponentially as 

Af{N) = Ar-7/2 gMoiV(i ^ 0{1/N^)), (93) 
where /Iq is known. The sum (4) can now be written as 

F{a,, a2, (3) = J2 e-^*^*-^*+<'^*+" {Nt+a\f\ N^) ^Af{Nt)Af{Nt+a) ■ (94) 

Nt,Nt+a 



The exponential part of ^^/Af{Nt)Af{Nt_^.a) can be absorbed in additive renormaliza- 
tions of the boundary cosmological constants Zt and Zt+a (i-e. additive renormaliza- 
tions of logctj, recall (4)). It follows that in the scahng hmit, i.e. for large N where 
we can use (93), the Laplace transforms of the matrix elements (A^i|T|A^2) a-re equal 
to the "7/2" fractional derivative'' of the free energy F{ai, a2, /?), 

(o o \ 7/4 
Q-g^j F{a,,a2,P). (95) 

where = jpA — as is clear from (5). 

oz a log Of ^ ' 

The scaling limit, and thus the continuum physics, is determined by the singular 
part of the free energy The leading behaviour of this singular part when wc approach 

^There are standard ways to define the concept of a fractional derivative, see for instance [13]. 
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a critical point as described in the previous section, is given by F{ai,a2, P) oc 
(5q;)^/^. It is now straightforward to apply (95) and one finds 



Q \ 7/4 / Q N 7/4 ^ 

-) F(q;i,q;2,/9) ~ (96) 
2 J oa 



d log aij \d log a2 

This is exactly the leading-order behaviour we expect for the transfer matrix when 
a — > from (1): 

{N^\e-''"\N2) ^ {N^\i\N2) = 5m,,n„ (97) 
and the Laplace transform of Sn-^^^n^ is (for large A'"'s) 

^ ^ :i + 0(a)). (98) 



5q;i + 5q;2 a^{Zi + Z2) 

It would be very interesting to expand to next order in a and thus obtain infor- 
mation about H. These terms come from the 0(a'^)-terms discussed in the previous 
section. This will only give us information about the matrix elements related to the 
states of the form (91), but as discussed in detail in [1] we expect that this is the 
only information relevant in the continuum limit of 3D gravity if the topology of the 
spatial shces is that of a sphere. 



Acknowledgement 

All authors acknowledge support by the EU network on "Discrete Random Geom- 
etry", grant HPRN-CT-1999-00161. In addition, J. A. and J.J. were supported by 
"MaPhySto" , the Center of Mathematical Physics and Stochastics, financed by the 
National Danish Research Foundation. J.J. acknowledges support by Pohsh Com- 
mittee for Scientific Research (KBN) grant 2 P03B 096 22 (2002-2004). 



References 

[1] J. Ambj0rn, J. Jurkiewicz, R. Loll and G. Vernizzi: Lorentzian 3D gravity with 
wormholes via matrix models, JHEP 0109 (2001) 022 [hep-th/0106082]. 

[2] J. Ambj0rn, J. Jurkiewicz and R. Loll: A nonperturbative Lorentzian path in- 
tegral for gravity, Phys. Rev. Lett. 85 (2000) 924-927 [hep-th/0002050]. 

[3] J. Ambj0rn, J. Jurkiewicz and R. Loll: Dynamically triangulating Lorentzian 
quantum gravity, Nucl. Phys. B 610 (2001) 347-382 [hep-th/0105267]. 

[4] J. Ambj0rn and R. Loll: Non-perturbative Lorentzian quantum gravity, causality 
and topology change, Nucl. Phys. B 536 (1998) 407-434 [hep-th/9805108]. 



20 



[5] J. Ambj0rn, J. Jurkiewicz and R. Loll: Nonperturbative 3-d lorentzian quantum 
gravity, Phys. Rev. D 64 (2001) 044011 [hep-th/0011276]. 

[6] J. Ambj0rn, J. Jurkiewicz and R. Loll: Computer simulations of 3d Lorentzian 
quantum gravity, Nucl. Phys. Proc. Suppl. 94 (2001) 689-692 [hep-lat/0011055]; 
3d Lorentzian, dynamically triangulated quantum gravity, Nucl. Phys. Proc. 
Suppl. 106 (2002) 980-982 [hcp-lat/0201013]. 

[7] V.A. Kazakov and P. Zinn- Justin: Two matrix model with ABAB interaction., 
Nucl. Phys. B 546 (1999) 647-668 [hep-th/9808043]. 

[8] V.A. Kazakov, M. Staudacher and T. Wyntcr: Exact solution of discrete two- 
dimensional B? gravity, Nucl. Phys. B 471, (1996) 309-333 [hep-th/9601069]; 
Almost flat planar diagrams, Commun. Math. Phys. 179 (1996) 235-256 [hep- 
th/9506174]; Character expansion methods for matrix models of dually weighted 
graphs, Commun. Math. Phys. 177 (1996) 451-468 [hep-th/9502132]; P. Zinn- 
Justin: Random hermitian matrices in an external field, Nucl. Phys. B 497 
(1997) 725-732 [cond-mat/9703033]. 

[9] J. Ambj0rn, J. Jurkiewicz and R. Loll: Renormalization of 3d quantum gravity 
from matrix models, [hep-th/0307263]. 

[10] P. Zinn- Justin: The asymmetric ABAB matrix model [hep-th/0308132]. 

[11] J. Ambj0rn, J. C orreia, C. Kristjansen and R. Loll: The relation between Eu- 
clidean and Lorentzian 2D quantum gravity, Phys. Lett. B 475 (2000) 24-32 
[hep-th/9912267]. 

[12] N.I. Muskhelishvili: Singular Integral Equations, Dover Publn. 1992. 

[13] S. G. Samko, A. A. Kilbas and O. I. Marichev: Fractional integrals and deriva- 
tives: theory and applications, Gordon & Breach, 1987. 

[14] S. Carlip: Quantum gravity in 2-hl dimensions, Cambridge Univ. Press, Cam- 
bridge, UK (1998). 



21 



